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Abstract 

We study the hydrodynamics of a freely-standing smectic-A film in the 
isothermal, incompressible limit theoretically by analyzing the linearized hy- 
drodynamic equations of motion with proper boundary conditions. The dy- 
namic properties for the system can be obtained from the response functions 
for the free surfaces. Permeation is included and its importance near the free 
surfaces is discussed. The hydrodynamic mode structure for the dynamics of 
the system is compared with that of bulk systems. We show that to describe 
the dynamic correlation functions for the system, in general, it is necessary 
to consider the smectic layer displacement u and the velocity normal to the 
layers, v z , together. Finally, our analysis also provides a basis for the theoret- 
ical study of the off-equilibrium dynamics of freely-standing smectic-A films. 
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I. INTRODUCTION 



A three dimensional smectic-A phase is at its lower critical dimension: the smectic 
layer displacement fluctuations grow logarithmically with the size of the system, and the 
divergence is sufficiently weak that finite size effects stablize laboratory samples 0. 

Smectic liquid crystals also form freely standing films with a controlled size and high 
degree of uniformity. Due to the presence of two free surfaces and finite thickness, ex- 
perimental and and theoretical studies on the equilibrium , and off-equilibrium |3| 
properties of freely standing smectic-A films show drastically different properties from bulk 
systems. Hence freely standing smectic-A films are excellent materials to study the effect of 
reduced dimensionality and surface effects. 

In addition to the static properties, the dynamics of a freely standing smectic-A film 
have been studied in the "low" frequency (u <C 50kHz) and "long" wavelength (A ^> 0.1 
cm) regime, where it behaves like a simple elastic membrane 0. The dynamic correlation 
functions have also been considered in the strongly overdamped regime, || both theoretically 
and experimentally. Notice, however, that these experiments were performed in a region 
where the effects of bulk elasticity are negligible; further experiments are needed to probe 
the dynamical crossover behavior of the system. 

In this article we introduce a continuum theory for the dynamics of a smectic-A film in 
the isothermal, incompressible limit based on the linear hydrodynamic theory constructed 
by Martin et al ||. Our study provides a finite-thickness counterpart to the hydrodynamic 
theory for bulk smectic-A systems. The following two questions are addressed and answered 
in our theory: (i) Permeation plays a role near the free surfaces ||, but how does it modify 
the smectic layer dynamics near the boundaries? (ii) The smectic layer displacement u, 
and the velocity of the liquid crystal normal to the layers, v z , are coupled hydrodynamic 
variables, but how do the hydrodynamic modes for u and v z for a film compare to the linear 
theory for a bulk smectic-A ? In the following sections we will show that the dynamic 
properties of the smectic layers in a freely standing film can be extracted from the response 
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functions of the free surfaces. Since the layer displacement for a smectic-A film has to satisfy 
certain boundary conditions on the free surfaces, the normal modes for the dynamics of the 
layers depend on both surface and bulk properties. We find that among all the normal 
modes, two are most strongly affected by the surface tension and film thickness. The decay 
rates for all other normal modes show the interplay of bulk-like behavior and effects from 
the surfaces. When the thickness of the film is small, the dominant thermal fluctuations of 
the layers come from these two modes, and the system can behave similarly to an ordinary 
fluid film. When the thickness becomes large, the other normal modes have sufficiently 
large contributions to the layer fluctuations, and the internal layer structure of the system 
can be detected experimentally. The present analysis reveals how the dynamics of a finite 
thickness smectic-A system depends on surface tension, viscosity, and also bulk elasticity. 
We also provide a theoretical picture for the crossover behavior of a freely standing smectic- 
A film from a thin film, which depends only on the surface properties, to a three-dimensional 
system with layer structure. Furthermore, this theory is also a starting point for theoretical 
calculations for the off-equilibrium dynamics of a freely standing smectic-A film. 

The remainder of this paper is organized as follows. Section II contains a review of the 
linear elastic and hydrodynamic theory for a bulk smectic-A in the isothermal, incompressible 
limit. After a brief review of the linear response functions of the surfaces of a smectic-A 
film, we extract the smectic layer dynamics and compare the bulk and film hydrodynamics in 
Section III. The crossover of smectic dynamics from film-like behavior to bulk-like behavior 
is discussed in Section IV. In Section V we summarize results and add concluding remarks. 
Some additional details on permeation near the surfaces and the calculation of the dynamic 
correlation functions are provided in Appendices A and B, respectively. 

II. ELASTICITY AND HYDRODYNAMICS OF BULK SMECTIC-A PHASE 

The symmetry broken direction for a smectic-A system is chosen to be the z-direction. 
Assume that temperature and density fluctuations occur on time scales much faster than the 
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characteristic times under consideration. In this limit we have an isothermal, incompressible 
system, and the elastic free energy of the smectic layer distortion for a bulk smectic-A in 
thermal equilibrium is provided to leading order by Jl[ 

1 du \ 2 „ / d 2 u d 2 u s 2 



where u(r, t) is the layer displacement, r = (x, y, z), and B and Ki are, respectively, the 
layer compression and undulation elastic moduli. The volume integral is performed over the 
bulk. In this article we consider exclusively a thermotropic smectic-A, for which typically 



B ~ 2.5 x 10 7 dyn/cm 2 , K\ ~ 1CT 6 dyn. The characteristic length A = JK\jB is typically 
on the order of the layer spacing (~ 10~ 7 cm) H. It has the following physical interpretation: 
in bulk systems the dominant thermal fluctuations of the layers with momentum q are those 
with q z ~ where q± is the magnitude of the projection of q on the xy-p\a.ne. When 
q± <^ A -1 , the dominant layer fluctuations come from q z <^ q±. 

The dynamics of a bulk smectic-A have been studied thoroughly in, for example, Ref . |§ , 
which we follow. In the absence of topological defects, under constant temperature and 
assuming incompressibility, to linear order the system satisfies the equations of motion 



Ov ■ 

p-Qj- = -diP + dja-j + h 5 iz (2) 



and 

du 

— = v z + C P h , (3) 

where p is the density (typically p ~ 1 g/cm 3 ), Vi is the ith component of the velocity field, 
and the pressure, p, is actually a Lagrange multiplier for the incompressibility condition. The 
viscous stress tensor is denoted by a' ||, ( p is the permeation constant, and the molecular 
field h is defined by 

»-*(£)■ 

In Eqs. (0) and (f|) we sum on repeated indices, and dj = d/dxj. For layer displacements of 
the form 



u(r, t) = u(q, u )j(<u--*->-+9'*)-** (5) 
with q z q±, permeation may be neglected [|]|§, and the relation between u and q is given 

by 



iL ° = ^ X i 1 ± ^-^M + ^l)) = iw±(ftL, e.) , (6) 

where 773 (typically 773 ~ lp) is one of the five viscosity coefficients of a smectic-A |§. 
The fact that u(r,t) couples to v z (r,t) in the equations of motion means that, in general, 
a hydrodynamic normal mode described by Eq. (Q is a linear combination of these two 
variables. There are two other normal modes for a given q, but they correspond to the 
diffusive modes of v x and v y , and they are decoupled from the motion of the smectic lay- 
ers [[H|. We will not discuss these two other modes in the rest of this section. In the limit 
4p(Bq 2 + Kiq\) / rj 2 q]_ <C 1, the viscous loss of the system is large, and the motion of the 
layers is strongly overdamped. Then iu is real, and the two decay rates are given by 

j_ = mol 
n p 

\_ = Bq 2 + K ig j 

T 2 vsql 

where 1/ti ^> l/r 2 . We can identify 1/ti as the decay rate of the diffusive mode associated 
with v z and 1/V2 as the decay rate for the smectic layer undulations. Since we are interested 
in the low frequency behavior of the system, it is appropriate to neglect the effect of T\. In 
the opposite limit, Ap(Bq 2 z + Kiq\_)/rj 2 q\ ^> 1, the elastic restoring force is much stronger 
than the viscous loss, and the system performs underdamped oscillations, i.e., iu becomes 
complex. The decay rate is simply t]^q'j_/2p, and the oscillation frequency is approximately 



(Bq 2 z + Kiqj_) / p , which is the frequency of "second sound" |l| . In this case, then, there are 
two propagating hydrodynamic modes associated with each q. Notice that when q z ~ Ag^, 
the system is strongly damped as long as pKi/rj 2 <C 1; for typical materials pKi/rj 2 ~ 1CT 6 , 
so that the dominant thermal fluctuations of the layers are overdamped. 
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III. HYDRODYNAMIC NORMAL MODES OF A FREELY STANDING 

SMECTIC-A FILM 



Let us now turn to the case of freely standing smectic-A films. In the absence of external 
fields, the smectic layers are always aligned and parallel to both free surfaces. We consider 
a film which extends from z = d/2 to z = —d/2 in the vertical direction and is of infinite 
extent in the x and ^/-directions. The geometry of the system is shown in Fig. 1, which shows 
as the displacement of the upper (lower) surface from its equilibrium position. In this 
section we first briefly review previous results |[TI] on the derivation of the linear response 



function for the equilibrium surface fluctuations of a freely standing smectic-A film in the 



low frequency, long wavelength regime [|T2|jL3|1 , then extend these results to the structure of 



hydrodynamic normal modes of the system. 

Suppose the system were perturbed by external forces with fixed in-plane wave vector 
and frequency u on both upper and lower surfaces, i.e, 

P e + Xt (r±,t) = P e + Xt (^uj) e^' r — * 

P- Xt (r±, t) = P- t (q±, lj) e^' r — ' . (8) 

We look for the linear response functions of the surfaces in the regime of weak external 
forces. On the free surfaces the system has to satisfy the following boundary conditions: 

be 1 The velocities of the free surfaces are the same as the normal component of the liquid 
crystal velocity on the surfaces. 



be 2 For free surfaces the normal component of the permeation force should vanish. |TT 
be 3 The force acting on the system is continuous across the free surfaces. 
The boundary conditions lead to the following equations for the system, 

—Qf~ = [««U/2 • (9) 
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du 

dz 



z=±d/2 



[diV z + d z Vi} z=±d/2 = 



dv 7 



du 



± _ p+2))3 _ +B _ 



2 a± 



«vic 



p 



ext 



i = x, y 
= 

z=±d/2 



(10) 

(11) 

(12) 



where a is the surface tension (typically a ~ 30 dyn/cm). 

Solving the equations of motion, i.e., Eqs. (0), and (^|), with these boundary conditions, 
one finds that, in the presence of the driving forces with given q_|_ and lu, there are four 
different q z s. One of them is associated with the diffusive motion for v x and v y , and is 
decoupled from the motion of u and v z ; hence it is irrelevant for the calculation of surface 
response functions (see be 1). Another two values of q z are associated with permeation, 
and these contributions to the layer displacement are small (see Appendix A). The last of 
these g 2 -values dominates the contribution to the smectic layer displacement; hence u can 
be approximated by 



u(r,t) = u s (q±,uj) cos [q z (q±,u)z] x 
+u A (q ± ,uj) sin[q z (q ± ,uj)z} x 



(13) 



Using ui±(q±, q z ) as defined in Eq.fl^), q z (q±,u) is determined by the solution of the following 
with Reg 2 > 0: 



iu = iu±(qj_,q z 



(14) 



where \q z \ <C qj_, and, in general, q z is complex under the conditions considered Recall 
that to and q± describe the driving forces and q z is not independent. This differs from the 
interpretation of the previous section in which the bulk hydrodynamics of a smectic-A was 
considered in the absence of external driving forces. 

When P^ xt = P~ xt {Pext = ~Pext)i t ne dynamics of the system are symmetric (antisym- 
metric) under z — > —z, and the layer displacement is described by us (ua) alone. We define 
the response functions for both symmetric and anti-symmetric surface motion as 



( S {q±, w) = X s (q ± , w) [P+t{q±, w) + P ext {q±, w, 

( A (q ± , w) = X A (q ± , w) \P e + xt (q±, w) - P- t (q x , w, ) , (15) 



where 



C S{A) = \ (C + ± C) ■ (16) 

A calculation from the continuum hydrodynamics leads to the following expressions for the 



response functions |TT 



X s = - - (17) 

2 aq\ - B q z (q ± ,uj) tan (rf q z (q±,u)/2) 

X A = I \ 

2 aq\ + B q z (q ± ,uj) cot (d q z (q±,u)/2) ' ' 

It appears that the permeation processes have no contribution to dynamics of the surfaces 
in the regime [|l2j where our analysis is done. However, the solution for u(r,t) in Eq. (|13D 



does not strictly satisfy be 2. This means that the contribution from permeation enables 
the system to satisfy the boundary conditions 0] but otherwise has little significance. The 
role of permeation in the dynamics of the system is discussed in further detail in Appendix 
A. 

We have derived the response functions X s< - A ' ) as functions of in-plane wave vector q±_ 
and frequency uo of the applied external forces. We can also define u on the complex plane; 
then the poles of the response functions provide the frequencies of the surface fluctuations 
in the absence of driving forces ||10|| . In linear theory these poles also reveal the frequencies 
of layer displacement. Each pole in the complex u;-plane corresponds to a hydrodynamic 
normal mode associated with the layer displacement. Since we are not interested in the 
dynamics of v x , and v y , which decouple from the dynamics of u and v z , these poles provide 
the information of interest. Two classes of normal modes with opposite symmetry properties 
can be found. One class of normal modes is symmetric under z —z; the frequencies are 
the roots of, 

Bq ^ u) = cot . (19) 



The other class of normal modes is anti-symmetric under z — > —z, and their frequencies are 
determined by 

M^l = _ tan t o.faMA . (20) 

aq{ \ 2 J 

Actually, we can go back and treat the above two equations as equations for q z ; we can find 
the z dependence of the layer displacement for the normal modes from these equations. The 
form of the layer displacement in Eq. (|l3f) suggests that, for given qj_, the layer displacement 
for an normal mode is an eigenfunction of the differential equation 

d 2 tp(q ± ,z 



dz 

with boundary condition 
bc2' 



-0(q±,*) (21) 



B — — ±aq ± V(q±, z) 



= , (22) 

z=±d/2 



It is straightforward to show that these eigenf unctions form a complete basis for the interval 
— d/2 < z < d/2. The apparent contradiction between be 2' and be 2 is discussed in 
Appendix A. 

In general q z in Eqs. ([T^) and ( p0| ) cannot be determined analytically. However, the q z s 
are all real, and we can choose all of them to be positive. The graphical solutions for the 
q z 's are shown in Fig. 2. We label the q z s for the normal modes as q z n \q±), n = 0, 1, 2, .... 
with q z °\q±) < q z l \q±) < q z 2 \qi_) < and note there are two normal modes associated 
with each q^ n \ The layer displacement for normal modes is even (odd) under z —z if n is 
even (odd). The layer displacement for the normal modes with n = is even under z — > —z. 
When both modes are underdamped, they form a pair of propagating waves, which in bulk 
systems corresponds to second sound |j] . When they are overdamped, they are two diffusive 
modes. In the strongly overdamped limit, the characteristic time scales are well separated; 
the fast mode corresponds to the diffusive motion of v z , and the slow mode corresponds to 
the layer undulation mode in bulk systems. 

To compare with the behavior of a bulk smectic-A, recall that translational invariance 
of the bulk system allows one to label modes by a (^-dimensional wave vector. For a freely 
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standing smectic-A film, the boundary conditions and the equations of motion select a set 
of normal modes for the dynamics of the system. We can express the normalized time- 
independent part of the normal modes in the following way, 

^(q_L, n; z) = N s (q±, n) cos [g^ n) (g ± )z] , (23) 

for n = 0, 2, 4, and 

V>(q±, n; z) = N A (q ± , n) sin \qi n \q ± )z\ , (24) 



for n = 1, 3, 5, ... , where 



[2 / sm(q^d)\ 1/2 



These normal modes satisfy the orthonormality condition 

rd/2 

/ dz ip(q_ L} n;z)4>(cL U m]z) = 5 mn . (26) 

J-d/2 

Hence, these modes can be used as a basis to expand any u(q±, z, t) as, 

«(q±, z, t) = X] w (q.u n > *) x ^(<u, ™; • ( 27 ) 

n 

From these relations one finds the dynamic correlation function for u. As discussed in 
Appendix B, the general form for the dynamic correlation function is 

c(q±,n,t) 

= < u(q±,n,t) u{— q±,m, 0) > 

iu+e-*"-* - luo-e-^ 1 k B T 

— Onm : : x TTo > l^°J 

where u;± = cu±(q±, q^) is related to q± and g( n ) through Eq. (^). In the underdamped case, 
u;± are complex and u;_ = uj* + . The dynamic correlation function has a damping rate r] 3 q"j_/2p 
and oscillates with frequency J —(^Qi/^P) 2 + {Bqi^ + i^igjj / p. In the overdamped case, 
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two decay rates exist, and the fast mode can be neglected only when the time scales are 
well separated. In general, when an experiment is performed outside the strongly over- 
damped regime, two characteristic time scales should be detected. For weakly overdamped 
motion, these are the two decay rates, while for underdamped motion, these are the real 
and imaginary parts of ou±. 

In summary, we find that in common with analysis of a liquid film, the normal modes are 
symmetric/anti-symmetric under z — > —z. They also are labeled by (q_i_,n), and the layer 
displacement can be treated as eigenfunctions of a differential equation. The frequencies for 
those normal modes can be obtained from the poles of the surface response functions in the 
complex oplane. As will become clear in the next section, the difference between those q^> 's 
which correspond to the normal modes for a finite thickness smectic-A film, and the set of 
numbers nn/d, n = 1,2,3, ... provides a measure of the degree to which finite thickness and 
surface tension change the dynamic properties of the layer displacement. Experimentally, 
these effects are typically obtained by measuring the frequencies and the magnitude of the 
autocorrelation function of the layer displacement. Layer dynamics well be discussed further 
in the following section. 

IV. CROSSOVER OF LAYER DYNAMICS FOR SMECTIC-A FILMS 

There are two special regimes where the q^'s can be evaluated approximately, and the 
physics of the system is also of interest. We call these the thin film and thick film regimes, 
respectively. We will show that the pair of normal modes with qf> is special since they 
strongly depend on the surface properties and show drastically different behavior from the 
thin film to thick film regime. 

When the dimensionless measure of film thickness 5 = aq\d/B <C 1, the solution for q z 
can be approximated by 

rm V2S [2a 

& = — = im (29) 

and 
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^ = t( 1 + w) = (t + ^IJ ■ "=1.2.3,- PO) 

For n > the modes are spaced approximately by n/d. This is the thin film regime (5 < 1). 

The equilibrium thermal fluctuations of the smectic layers in this limit are dominated by 
surface effects. This can be illustrated by considering the pair correlation function for the 
layer displacement. From Appendix A, in thermal equilibrium, for a smectic-A film with 
given q±, the contribution to the pair correlaton function from the modes with is 

(K<U,0)| 2 )= *f 4 . (31) 

Furthermore, when K\q\dja <C 1, the contribution from layer elasticity is negligible, and 
the thermal fluctuations of the smectic layers for these normal modes are independent of the 
bulk elastic coefficients ||. For a typical material, B ~ 2.5 x 10 7 dyn/cm 2 , K\ ~ 10 -6 dyn, 
a ~ 30 dyn/cm, and the thin film regime corresponds to q\d < 5/« ~ 10 6 cm -1 . The layer 
fluctuation is dominated by the surface tension term when q\d <C ol/K\ ~ 3 x 10 7 cm -1 , 
i.e., in the thin film regime, layer fluctuations due to these two normal modes are typically 
dominated by the surface tension term. In the following we will simply drop the contribution 
from the K\ elasticity term in the thin film regime. 

The total contribution from the fluctuations of the other normal modes for fixed q± is 
approximately 

j2 ^Z. 

n= i B (nn/d + 2aq\/mrB) 
~ -f - ■ (32) 
This contribution is small compared to the contribution from the normal modes with q^ in 
thin film limit. As a result, we find that a freely standing smectic-A film can be treated as 
an ordinary fluid film with surface tension a and thickness d in this limit. || The internal 
structure of the layers is not important. Comparing with bulk smectic-A systems where the 
dominant layer fluctuations for given q± come from q z ~ \q\ , for a thin film the dominant 
layer fluctuations come from only two normal modes, and it is in general not true that 

# ~ Ml- 
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From Eq. @ the normal mode frequencies u)±(q±, q^) are given by 



\ 2 



l-^jl, (33) 

where g c = \JSpa/rjld is a characteristic wave number of the film. In the case of short 
wavelength fluctuations with q± ^> q c , viscous loss is strong compared to elastic effects, and 
the overdamped decay rates are 

■ (o)/ \ VsQl 

^V) = | . (34, 

Since we are interested in slow motion of the system, iuj_ (q±) provides the characteristic 
time for the system. Previous experimental studies of dynamic correlation functions were 
performed in this regime || . Notice that in this regime the decay rate is independent of the 
magnitude of the in-plane wave vector and bulk elastic constants. In general, experiments 
which probe the fluctuations of the film with smaller in-plane wave vector should be able 
to identify two decay rates. This is especially so when the motion of the system is close to 
being critically damped, as can be seen from Fig. 3. For long wavelength q^ <C q c the film 



oscillates with frequency y2a/pd q±_ and damping rate r] 3 qj_/2p. Again, both time scales in 
this range of q± are independent of B and K\. As q± — > 0, the damping is small; this is the 
regime in which the long wavelength theory and experiments |7j have been performed. 

We have shown that the physical properties for the system in the thin film regime are 
basically the same as a two dimensional object, i.e., a simple fluid film. The dynamics of 
the system show a competition between the elasticity due to surface tension and viscous loss 
characterized by the coefficient 773. The underdamped motion for a smectic-A film should 
not be confused with "second sound" for smectic layers, because for the thin film the layer 
compression is negligible. In this regime, the underdamped motion is simply the vibration 
of a simple fluid film with surface tension a. Also notice that although the magnitude of 
q(°^ depends on the elastic constant B, this dependence does not show up in the dynamical 
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correlation function for the hydro dynamic normal modes, as can be seen in Eq. (|3T|), which 
is the same as that for the surface fluctuation. 

To probe the crossover behavior of the film, experiments have to be performed outside the 
thin film regime so that the effect of the bulk elasticity can be revealed. This can be done by 
increasing the film thickness or adjusting the in-plane wave vector. When aq\d/B increases 
to order unity, the layer structure of the system starts to play a non-negligible role in the 
system. First of all, this is reflected in the fact that the approximation for q^ in Eq. (gD 
is poor, and the higher order terms contribute. Then the n = modes start to show the 
influence of bulk elasticity. Second, the layer fluctuations of all other normal modes become 
comparable to the layer fluctuations for the n = modes. Hence, the dynamics of the 
layers for fixed q± are no longer simply controlled by u + (q±, q^) and tu_(q±, q^), and this 
means that there are more relevant degrees of freedom in the system than the displacement 
of the two surfaces. The first sign of this crossover behavior can be detected by exciting 
the n = 1 modes. As pointed out in the study of the surface dynamics of freely standing 
smectic-A films |TT|jr4| , this pair of modes is similar to the so-called peristaltic mode for a 



soap film. However, the characteristic time scales of these modes strongly depend on the 
elasticity of the smectic layers; hence, the dynamic properties of these modes actually reveal 
the difference between a smectic film and a soap film. Possible future surface light scattering 



experiments for probing these normal modes are also discussed in Ref. [TT[ . 

For 5 1, qf^ increases and q^d ~ 7r; the q^s for all normal modes can be approximated 
by a single expression, 

(n) (n + l)ir / 2\ (ra + lW/ 2B \ 

From Eqs. ( |23"D and ( |24] ) we find that the layer displacement for the normal modes on both 
free surfaces becomes vanishingly small as 5 _1 — > 0. This means that the layer displacement 
is not affected by the surface tension in this limit, and the behavior of the system is the 
same as in a bulk system. Hence, the thick film regime is defined by 5 = aq]_d/ B 3> 1 
The dynamics of the layers in the thick film regime is similar to that of a bulk smectic-A. 
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For fixed q±, when q^> is large such that Ap(B q^> +K 1 qj}/r)^q'^ > 1, the smectic layers are 
underdamped. This is basically the second sound of the layers, but with a correction due to 
surface tension and finite film thickness. When 4p(_B q^ 2 + Kiq\) /rfeq\_ <C 1 the motion of 
the layers is strongly overdamped. The decay rates which are slow and correspond to layer 
undulation are 

. (n) Bq^ + K.qj 

iw_(q ± ,q { z >) = 2 > ( 36 ) 

V3Q± 

while the fast decay modes now correspond to the diffusive mode for v z . 

To study the dynamics of the smectic-A film in the thick film regime for typical laboratory 
materials, (B ~ 2.5 x 10 7 dyn/cm 2 , a ~ 30 dyn/cm) we need q\d 3> B/a ~ 10 6 cm -1 
(5 ^> 1). For d ~ 10~ 3 cm, which can be achieved easily, the momentum transfer in the 
xy-plane should be at least of the order 10 5 cm -1 , which can also be achieved. 

The crossover of the dynamics for the system is illustrated in Fig. 4, where the normalized 
dynamic correlation functions C(q±,n,t) = c(q±, n, t)/c(q±, 0, 0) for normal modes with 
n =0, 1, 2, 3, are shown for a smectic-A film with thickness d = 10 /iin, and typical elastic 
coefficients and viscosities. The in-plane wave vector q± is chosen such that the dimensionless 
parameter = Bjaq\d is 20000, 2000,1, and 0.1, respectively. When = 20000, the 
normal modes with n > are negligible and n = normal modes are underdamped. For 
5 1 = 2000 the n > normal modes are still negligible but the n = normal modes are 
overdamped with a decay rate which is not a single exponential. This occurs when the two 
time scales are close to each other, as illustrated in Fig. 3. For 5 = 1 the normal modes 
with n > can be observed; this is the crossover regime between the thin film and thick 
film limits. For <5 -1 = 0.1 the contribution from bulk elasticity is sufficiently large that the 
dynamic correlation functions for normal modes with n > are easily seen in the figure. 
Similar to the undulation mode in bulk systems, all normal modes decay exponentially. This 
is the overdamped limit of the thick film regime. The actual q± values in the figure change 
over many orders of magnitude. Thus, experimentally probing this entire figure requires 
different scattering probes. 
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V. CONCLUDING REMARKS 



The layer dynamics of a freely standing smectic-A film have been determined by analyzing 
the response functions for the surface displacement. By determining the positions of the 
poles of the response functions in the complex u-plane, we find not only the frequencies but 
also the spatial configurations for the hydrodynamic modes of the system. The fact that u 
and v z are two non-separable hydrodynamic variables, in general, is reflected in the form of 
the dynamic correlation function for u. 

When the dimensionless thickness 5 = aq\dj B -C 1, the internal structure of the layers 
is not important for the system. For fixed q± one pair of normal modes dominates the stat- 
ics and dynamics of the system. Furthermore, the layer fluctuations and the characteristic 
frequencies for these important normal modes show basically no dependence on the elas- 
ticity of the smectic layers, i.e., only surface tension and viscosity are important. To find 
information on the internal structure of the system, i.e., the smectic layers, it is necessary to 



excite other normal modes. As discussed in reference ||11|| , exciting the n = 1 modes provides 
information on the crossover behavior of the system from a quasi-two-dimensional system 
to a three-dimensional system. This can be done by adjusting the experimental parameters 
such that 5 = aq\djB is of order unity. The mode structure for the system in the regime 
5 ^> 1 is similar to that of bulk systems. However, the characteristic time scales for the 
hydrodynamic modes now acquire surface and finite thickness corrections as well. 

To summarize, we have explored the dynamical properties of a freely- standing smectic-A 
film in the linear regime. It is shown that existing experiments were performed in the limit 
where only the two-dimensional character of the film can be detected. Future experiments 
on the dynamical correlations of the smectic layers can reveal the crossover behavior of the 
system to a regime in which layer elasticity begins to play a role. Our work also provides 
the basis for a future theoretical study on the dynamics of a freely-standing smectic-A film 
far from equilibrium, where the system has been studied experimentally and shows behavior 
which is drastically different from that of a bulk system M . The effect of finite thickness and 
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surface tension will have to play a central role in a theory which describes such dynamical 
behavior. 



ACKNOWLEDGMENT 

We thank Professor X-l. Wu for helpful discussions. H-Y.C is grateful for fellowship 
support from the University of Pittsburgh. DJ is grateful for the support of the NSF under 
DMR9217935. 

APPENDIX A: 



We discuss permeation processes in a free-standing smectic-A film in this Appendix. We 
look for solution of the equations of motion with boundary conditions and external forces 
given in Section |I[ In the long wavelength, low frequency limit [12| the solution is of the 
form 



U = U ei + U p 



(Al) 



where 



u d (r,t) = u s (qj_,w) coa[q g (q±,u)z] x e iq± ' r± *"* 
+u A (q ± ,u) sm[q z (q x ,aj)z] x e ^ r ^«" 



(A2) 



is the approximate solution used in the text. It is independent of permeation constant. The 
second contribution 



u p (r,t) 



<(q±) x e +(w)p(i±)(z-d/2) +cc 



x e 



+ 



U, 



"(<U) X e -(l+0p(9±)(*+d/2) +CA 



x e 



(A3) 



contains an inverse length, p(q±) = y q±/2((pi]3) 1 / 2 , which characterizes the decay of per- 
meation away from the free surfaces. In typical materials p(q±) 3> q±_ S> q z in the low 
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frequency, long wavelength regime ||12|| . The amplitudes contain the contributions from 



permeation and are related to u e i through the full boundary conditions. 

From Eqs. ([P2]) and (|I0|) we find that \u^ | -C \u e i\, i.e., the contribution of permeation 
to the layer displacement is small ||11||. Then Eq. (0) leads to 



C ± « [u e l]z=±d/2 



(A4) 



It is straightforward to show that the elastic free energy is 



+ E 



dV { D 
1 



i=+, 



fdu\ 
[to) 

dS i a 



+ E 



dV < B 

dS l a 



i=+, 



foe; 

\ dx 
+ B 
( du,) 



' d 2 u d 2 u 
+ 



dx 2 dy 



dx 



' dug 

dz 



8y 

\ 2 



d 2 u e i d 2 u 



dx 2 



+ 



ti 



dy 2 



' du 



ti 



dy 



(A5) 



z=±d/2 

In the linear theory the contribution from permeation can be separated from the contribution 
from u e i. Since \u^ \ <C \u e i\ in the regime where we perform this long wavelength analysis, 
we can neglect the u p part. However, it is u e i, not the full displacement u, that can be 
expanded as a linear combination of the normal modes discussed in the text and satisfies 
be 2' (Eq. (|22"D ). When u v is included, the permeation force Bdu/dz vanishes on the free 
surfaces due to the existence of u p . This is noted in the literature ]15[ but has not been 
included in the previous study of the static and dynamic [|J theory for free standing 
smectic-A films. 



APPENDIX B: 

In this Appendix we calculate the dynamic correlation functions for u and v z . We consider 
the u e i part only; the u p part simply corresponds to the permeation mode discussed in the 
literature fl]||. Our starting point is the equations of motion, i.e., Eqs. (§) and (j^), and the 
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dynamics of the normal modes labeled by (q±,n) are considered. Eliminating the pressure 
via the incompressibility condition, we find that there are two diffusive modes associated 



with v x and v y , but v x and v y decouple from v z and u JT0|. Neglecting permeation, |TBj the 
equations for v z and u can be expressed as 



d_ 

Of 



where 



D 



1 \ 

-R -D 



mq 2 ± 



(Bl) 



(B2) 



represents the dissipative part, and 

R=- (Bq^ 2 + K iq \) (B3) 

represents the reactive part. We have used the fact that q^ <C q± for all n in the regime of 
our calculation to simplify the expressions. The eigenvectors of the above matrix equation 



are 



± 



1 



\ 



V 



-iu±{q±,qi n) ) 



(B4) 



with eigenvalues — iu±(q±, q^) given in Eq. @. Straightforward algebra leads to the fol- 
lowing relation 



v z (t) 



1 



V 



[iu + e — ito^e 



^+*]<u(0) + [e - * 6 '-* - e -<w + t ]u ig (0) 



(B5) 



tu + c<j_[e + — e 

where we used simplified notation, tu± = uj±(q±, q{ n ^), u(t) = u(q±,n,t), v z (t) = v z (q±,n,t), 
etc. From this one concludes that the dynamic correlation function for u is 

< u^Li^.t) u(—q±,m,0) > 

iuj+e-™-* - iu_e- iuJ+t . . , l2 
= Onm : : x < \u(q±,n, 0)| > 



ILOa 



iuo. 



x , (B6) 

Bq z n) +K iq \ 

as given in Section |TJ. The calculation for the correlation function for v z is similar to that 
for u. 
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Figure Captions 

Figure 1. Schematic of a freely standing smectic-A film of thickness d. The y— axis is 
pointing into the paper. The dotted lines are the equilibrium positions of the free surfaces. 

Figure 2. Graphical solution for the q^s for dimensionless parameter B/aq\d = 0.5. The 
abscissa in all cases is q z d/n. The three functions plotted are: Bq z /aq\ (dot-dashed line), 
cot(q z d/2) (solid lines), and —tan(q z d/2) (dashed lines). The crossings of the dot-dashed 
line and the solid lines give the values of is q^d/n for normal modes with n even; the 
crossings of the dot-dashed line and the dashed lines (except the origin) give the values of 
q^d/ir for normal modes with n odd. 

Figure 3. Decay rates for normal modes with n = in the thin film regime. The 
dimensionless wave vector is defined as Q = q^d. The material parameters are chosen as 
r] 3 = 1.0 p, p = 1.0 g/cm 3 , d — 10 /im, and a = 30 dyn/cm. When the decay rates merge, 
the system undergoes underdamped motion. 

Figure 4. Dimensionless dynamic correlation function C(q±, n, t) = c(q±, n, t)/c(q±, 0, 0), 
for n = 0,1, 2, 3. The material parameters are chosen to be d — 10 /im, B = 2.5 x 
10 7 dyn/cm 2 , K x = 10~ 6 dyn, a = 30 dyn/cm, p = 1.0 g/cm 3 , i] 3 = 1.0 poise, and q± is 
chosen such that S^ 1 = Bjaq\d = 20000, 2000, 1, and 0.1, respectively. 
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